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Abstract 

An analysis of one-loop vacuum fluctuations associated with a scalar fleld conflned in the 
interior of a inflnite waveguide of rectangular cross section is presented. We flrst consider the 
massless scalar fleld deflned in a four-dimensional Euclidean space. To identify the inflnities of the 
vacuum fluctuations we use a combination of dimensional and zeta function analytic regularization 
procedures. The inflnities which occur in the one-loop vacuum fluctuations fall into two distinct 
classes: local divergences that are renormalized by the introduction of the usual bulk counterterms. 



and surface and edge divergences that demand countertems concentrated on the boundaries. We 
present the detailed form of the surface and edge divergences. Finally we discuss how to generalize 
our calculations for a confined massive scalar field defined in a higher dimensional Euchdean space. 
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1 Introduction 

In 1948 Casimir predicted that uncharged, parallel, perfectly conducting plates should attract 
with a force per unit area, F{L) oc where L is the distance between the plates. This force 
can be interpreted as the manifestation of the zero-point energy of the electromagnetic field in the 
presence of the plates. Complete reviews of this effect can be found in refs. ||2| [^. 

The aim of this paper is to generalize the results obtained by Fulling |^ and others, for the 
renormalized vacuum fluctuations of a massless scalar field calculated between two parallel fiat 
mirrors. Our purpose is to analyse the vacuum fluctuations near surfaces where the field satisfies 
some classical boundary condition, the boundary being defined by two of the coordinates which we 
assume are closed intervals. A special case of this situation is the infinite waveguide of rectangular 
cross section. We are thus interested in obtaining the renormalized one-loop vacuum fluctuations 
associated with a scalar field defined in the interior of a infinite waveguide of rectangular cross 
section. We first consider a massless scalar field in a four-dimensional Euclidean space, and then 
extend our calculations to a massive scalar field defined in a higher dimensional Euclidean space. 
As stated by Deutsch and Candelas ^ and also Candelas ^ the calculation of the renormal- 
ized one-loop vacuum fluctuations contains all the important steps that we need to calculate the 
renormalized stress-energy tensor associated with the scalar field. 

It is well known that there are two quantities which might be expected to correspond to the 
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total renormalized energy of quantum fields H] . The first is the mode sum energy, defined as 



mode 



oo 

duj-uj[N{uj) - No{uj)], 



(1) 
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where is the zero-point energy for each mode, N{ijj)duj is the number of modes with frequencies 
between oj and uj + du m. the presence of boundaries, and NQ{uj)duj is the corresponding quantity 
evaluated in empty space. In this case the divergences that appear in the regularized energy 
are given by Weyl's theorem and its generalization, that relate the asymptoptic distribution of 
eigenvalues of some elliptic differential operator with geometric parameters of the surface where the 
fields satisfy some boundary condition 0. The second is the volume integral of the renormalized 



energy density (-E)™^ obtained by the Green's functions method [|T0| ||ri| ||T2[ fTS 

For special configurations where the modes and the eigenfrequencies of the electromagnetic 
field in the presence of the boundaries can be found, it is possible to find the Casimir energy and 
the Casimir force. For the case of cylindrical geometry, Balian and Duplantier and also others 



calculated the Casimir energy of the electromagnetic field [|T4| |[T5| [T^ [|l^]. For the important case 
of scalar or fermionic fields in the presence of a spherical shell. Bender and Hays |[13[, studying the 
global problem, found the renormalized zero-point energy of these fields assuming that the fields 
are confined in a spherical region of the space. Many years before Bender and Hays, Boyer [T^ and 



Davies [19 also studying the global problem, obtained the Casimir energy of an electromagnetic 
field in the presence of a perfectly conducting spherical shell. A systematic study of the spherical 



shell configuration was made by Milton He calculated the zero-point energy for gluons and 
fermions, assuming that the fields are confined in the interior of the shell. He later calculated 



the Casimir energy of massless fermions in the presence of a spherical shell, in this case taking 



into account the external modes of the fermions |21]. More recently, working in a generic flat 



d-dimensional spacetime Bender and Milton obtained the Casimir energy associated with a 
massless scalar field taking into account the interior and exterior contributions of the modes in the 



presence of the hypersphere. Still studying the spherical configuration, Romeo p3| investigated 
the Casimir energy of a massless scalar field and for QED assuming that the field is confined in 
the interior of the spherical shell. Finally Bordag et.al. have calculated the Casimir energy 
associated with a massive scalar field in the presence of a spherical shell assuming that the interior 
and the exterior modes give contributions to the energy. 

Although the main interest in the literature is the global approach, where the Casimir energy 
can be found and from which one can derive the force on the boundaries, the necessity of studying 
the local problem has often been suggested. First, the local properties of the vacuum fluctuations 
can in principle be observed by measuring the energy level shift of atoms interacting with an 



electromagnetic field |^. Second, the distortion of the vacuum fluctuations due to the 
presence of classical boundaries can also be measured by studying the spontaneous and induced 



emission of excited atoms in the presence of classical boundaries |2^. For a update discussion 
of QED vacuum effects, see for example ||2^. Finally, it is clear that local results contain more 
information than the global ones. A few years ago. Actor |]5D| and also Actor and Bender considered 



this kind of problem ||3T|. These authors studied the use of the generalized zeta function method 



2| 1 34 1 1 35] to find the one- loop effective action associated with a scalar field defined in 
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the interior of a infinite waveguide of rectangular cross section. The local problem has also been 
considered by other authors. For the case of parallel plates geometric configuration Brown and 
Maclay |jl2|, using the Green's function method obtained the renormalized vacuum expectation 
value of the stress-energy tensor associated with an electromagnetic field. Deutsch and Candelas ||^ 
evaluated the renormalized stress-energy tensors associated with a conformally coupled scalar field 
and also with an electromagnetic field in the wedge-shaped region formed by two plane boundaries. 
Recently, Brevik et al [^] repeated these calculations using Schwinger's source theory. Their 
results agree with those of Deutsch and Candelas. In the spherical geometry, the local problem 
was also investigated by Olaussen and Ravndal [^, who studied the vacuum fluctuations of 
an electromagnetic field within a perfectly conducting spherical cavity. They found that the 
vacuum expectation value of the squared electric and magnetic fields diverge as one approaches 
the boundary. This result has also been obtained by Deutsch and Candelas 0, DeWitt |3^, and 
Kennedy et al [^. Olaussen and Ravndal |^0| and also Milton generalized this result to the 
non-abelian gluon fields in the MIT bag model. It has often been suggested that these surface 
divergences are related to the uncertainly relation between the field and the canonical conjugate 



momentum associated with the field fS^] [@ @]. In other local calculations. Ford and Ford 



and Svaiter ||4^ discuss the possibility of amplification of the vacuum fluctuations. These authors 
studied the renormalized vacuum fluctuations associated with a scalar and electromagnetic field 
near the focus of a parabolic mirror. Using the geometric optics approximation they found that 
the parabolic mirror geometry can produce large vacuum fluctuations near the focus, much as 
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what happens in the classical focusing effect by the parabolic mirror geometry. 

In our study of the one-loop vacuum fluctuations associated with a scalar field we develop a 
method adequate to deal with rectangular geometries. In this paper, using analytic regularization 
procedures, we first calculate the regularized vacuum fluctuations associated with a massless scalar 
field, confined in the interior of a infinite waveguide with rectangular cross section, in a four- 
dimensional Euclidean space. We first rederive the well-known result that there are surface and 
edge divergences that require the introduction of surface |^ |^ and edge counterterms in 
the renormalization procedure. Then, we show how it is possible to generalize our results to 
a massive scalar field in a higher dimensional Euclidean space. Preliminary calculations of the 
renormalized vacuum expectation value of the stress-energy tensor in the rectangular waveguide 
were performed by Dowker and Banach [|^. Also {E)^°^^ in rectangular geometries has been 
calculated in refs. [0]. A seminal paper studying these geometries was made by Ambjorn 
and Wolfram [^. More recently Milton and Ng studied the Casimir effect in (2 + 1) Maxwell- 



Chern-Simons electrodynamics in a rectangular domain Also Hacyan et al ||5^ and Maclay 
]5l| respectively studied the vacuum fluctuations of the electromagnetic field and the Casimir force 
in the interior of a rectangular waveguide. 

When studying interacting field theories for translationaly invariant systems, one usually goes 
from coordinate to momentum space representation, a more convenient framework to analyse the 
divergences of the theory. In this representation translational invariance is expressed by momen- 
tum conservation conditions. Because the system of interest for this work possesses translational 
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invariance along two directions, it is more convenient to use a mixed coordinate-momentum rep- 
resentation for the Green's functions [E^ . For a recent treatment of systems without translational 



invariance, see for example and also As we discussed, the fundamental problem in the 
infinite waveguide of rectangular cross section is the lack of translational invariance, which man- 
ifests itself by the fact that the Green's functions associated with the scalar field are expressed 
in terms of infinite double summations. Although the one-loop vacuum fluctuation is written in 
terms of expressions involving double summations, a simple trigonometric identity allows us to 
obtain expressions with only one summation. The advantage is that all the calculations can then 
be done analytically. We show that the form of surface and edge counterterms that we have to 
introduce to renormalize the one-loop vacuum fluctuations, can be explicitly calculated and we 
discuss some different physical arguments that support the introduction of surface countertems 
that remove these divergences. 

Throughout the paper we use the term bulk counterterms to those conterterms that are not 
related to the existence of the boundaries. For instance, in the case of periodic boundary conditions 
we have these counterterms. As well as a second group that we call surface and edge counterterms 
that are related to the existence of the boundaries. 

The organization of the paper is the following: In the section II we study the one-loop vacuum 
fluctuations of a massive scalar field confined within a rectangular waveguide. In section III we 
use an analytic regularization method to identify the infinities that appear in the expression of 
the vacuum fiuctuations of the massless field near the boundaries. In section IV, we show how it 
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is possible to generalize our calculations to a massive scalar field defined in a higher dimensional 
Euclidean space. Conclusions are given in section V. Finally in the appendix we present the formal 
relation between the one-loop vacuum fluctuations and the local generalized zeta function. In this 
paper we use h = c = 1. 



2 Vacuum fluctuations of a scalar fleld conflned in a rect- 
angular waveguide 

In this section and in the next one we will investigate the one-loop vacuum fluctuations associated 
with a scalar field defined in the interior of a infinite waveguide with rectangular cross section. 
First, we assume that the scalar field is defined in a four-dimensional Euclidean space, where the 
last two coordinates are unbounded, while the first two, which we call Xi and X2, lie in the interval 
[0,a] and [0,6], respectively. We assume Dirichlet boundary conditions on the boundaries. The 
free field is defined in the region 

= X = {xi, X2, X3, X4,) : < xi < a, < X2 < b C R"^, (2) 

with Dirichlet boundary conditions at a;i = and xi = a, and also X2 = and X2 = b. As 
stressed previously, the lack of translational invariance introduces surfaces and edges divergences. 
One way to reduce the degree of these divergences is to smooth out the surface of the plates, 
for example by using soft, hard or semi-hard boundary conditions The one- loop 



graphs will depend on the ad-hoc model assumptions, and consequently we prefer to maintain the 



hard walls. Instead of smoothing the plates surfaces, another way to avoid surface divergences, 



discussed by Kennedy et al is to treat the boundary as a quantum mechanical object. This 
approach was used recently, by Ford and Svaiter in the case of parallel flat plates, to solve a long 
standing paradox concerning the renormalized energy of minimally and conformally coupled scalar 



fields As we stressed before, we prefer, at least at the moment to keep only a hard classical 
boundary conditions. 

We would like to stress that we are studying first the massless four-dimensional case since 
this simpler case will give us an indication of the behavior that may be expected in the most 
general case of a massive scalar field in a higher dimensional Euclidean space. A logical going 
is to use analytic regularization procedures to identify the divergent contributions that appear 
in the one-loop vacuum fiuctuations. Let us start using dimensional regularization, working at 
the begining in a d-dimensional Euclidean space. Assuming Dirichlet boundary conditions on the 
walls we calculate first the two-point Schwinger function at coincident points in the interior of 
the waveguide. Then we use analytic regularization procedures to identify the form of the surface 
and edge divergences. In the one-loop vacuum fiuctuations, {ip'^{xi,X2,a,b)), we will change the 
notation of the previous section to xi = x, X2 = y (note that in our notation {(p'^{x,y)) means 
the one-loop vacuum fiuctuations in one point with cartesian coordinates (x, y) of the rectangular 
cross section of the waveguide). Let the waveguide be oriented along the z axis with walls at x = 
and a and y = and b. In this mixed representation, since we are assuming Dirichlet boundary 
conditions, the expression for the one-loop vacuum fiuctuations, Tdd{x, y, a, b, d) = {(p'^{x, y)) can 
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be written as: 



TMx, y, a, b, d) = ^^^^ ^ sm (^) sm 



/ d'-'p-. (3) 



There are two points that we would hke to stress. First is the fact that to perform analytic 
regularizations, we have to introduce a parameter /i with dimension of mass in order to have 
dimensionless quantities raised to a complex power. For sake of simplicity, we omit the yU factor 
in the following. Second is the fact that the generalization for the case of Neumann boundary 
conditions is straightforward, although in this case infrared divergences associated with the n = 
mode will appear in the case of massless scalar field. To circumvent this situation we must have 
a finite Euclidean volume to regularize the model in the infrared. 

Using trigonometric identities, it is convenient to write the one-loop vacuum fluctuation in the 
following way : 

Tdd{x, y, a, b, d) = T{a, b, d) + T(x, a, b, d) + T(?/, a, 6, d) + T(x, a, 6, c?), (4) 

where each expression of the above equation are given by the following. For T(a, 6, d) we have 

4 °° /■ 1 
T{a,b,d) = ——r— y d'^-^p-, ^. (5) 

The second term in Eq.(^, T{x, a, b, d) is given by 

T(x,a,b,d) = —T{a,b,d) + ——r^y d'^-^p-, ^. (6) 
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The expression for T{y, a, b, d) has the same functional form of the above equation only changing 
X by y and a by h. Consequently we have : 

and finally : 



T{x,y,a,b,d) = — 
[Zn 



1 OO „ 1 

- y / c;rf-2p t 



( 2n7fx ,2n'ny 2mTX ,2n'7iy \ 
1 — cos( ) — COS[ — - — ) + cos( ) COS[ — - — ) I . (8j 

Thus the contribution given by Eq.(8) contains the same contributions given by Eqs.(5),(6),(7), 
as well as a contribution that contains edge divergences that we will define as iV(x, a, 6, d). It is 
given by 

N{x,y,a,b,d)= , , y d'^-^p^ , ' \ . (9) 

Let us study each contribution separately. Using dimensional regularization on Eq.(^) it is possible 
to write T(a, b, d) in terms of the Epstein zeta function. Thus we have 

^' ^) = (2 Jy-^ab ^^^ - ^) £ 7 ' , (10) 

{^^TT) aO / n,n'=l (^2 _^ (-"21)2 _^ fllljLy] 2 



a ' ^ b 

The contribution given by T(a, b, d) is one part of the vacuum field fluctuations which does not 
depend from the distance to the boundaries and in the renormalization procedure will require 
only a usual bulk counterterm. The form of the counterterm is given by the principal part of the 
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Laurent expansion of Eq. ( jlOD around some d, which must be given by the analytic extension of the 
Epstein zeta function in the complex d plane. The structure of the divergences of the Epstein zeta 



function is well know in the literature |6^. Since the polar structure of T(a, 6, d) can 



be found in the literature, to calculate the analytic structure of T{x, a, 6, d) we will concentrate 
only on the position dependent divergent part given by T(x, a, 6, d) + |r(a, 6, d). This expression 
is given by 

1 2 °° /■ cosf^^^) 

Although Eq.(pJ]) is written in terms of two sums, one of the sums can be easily done using a 
trigonometric expression given by p3 p3 : 



^ cosnt 1 TT coshy4(7r — t) 

}^^n'^ + A^^~2^^2A sinhTTTl ' ^^^^ 



Using Eq. ([T2|) , it is possible to write Eq.([TT|) as 



where : 



T{x, a, 6, d) + 2-^('^' ^' ~ Ri{0', b, d) + R2{x, a, 6, c?) (13) 



1 oo „ 1 



and 



1 ^ r , , 1 cosh((a-2x)Jp2 + m2 + (^^)2) 



(27r)" " 6 + m? + {^f sinh(a Jp^ + ^2 + (n^\2^ 
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It is clear that to calculate the analytic structure for the case of the position dependent divergent 
part T{y,a,b,d) we can use the same method that we use for T{x,a,b,d). Consequently the 
expression for T{y, a, b, d) + |T(a, 6, d) is : 

T{y, a, 6, d) + ]jT{a, b, d) = /i(a, b, d) + ^{y, a, 6, d) (16) 

where : 
and 



l2{y,a,b,d) = / / d p—^ , . (18) 

{2n)'^-\^J + m2 + (^)2 sinh(6Jp2 + ^2 + (1^)2) 



Using dimensional regularization in both expressions /i(a, 6, d) and /2i(a, 6, (i), we obtain a special 
Epstein- Hurwitz zeta function. The analytic extension of this function for general d in the massive 
and massless case can be found in the literature. For the massive case see, for example ||6^. For 
the massless case, the duplication formula for the Gamma function allows us to write 

/i(a, b, d)\m=o = -^h[d)a^ - d)V{A - d), (19) 



where 

is an entire function of d. Since ({z) can be analytically continued from an open connected set of 
points in the complex z plane into the entire domain of z, it is easy to find the analytic structure 
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of Ii{a,b,d)\m=o in the complex d plane. It is clear that Ii{a,b,d) and Ri{a,b,d) demand bulk 
counterterms. To find the analytic structure of hiu, a, b, d) and R2{x, a, 6, d), let us concentrate on 
l2{y,a,b,d). Integrating over the solid angle in Eg. (pISD , i.e using the fact that d'^p = p'^'^dpdQd 

d 

and / dVLd = Itzt we have : 



hiy,a,b,d) = -hid)y^ dpp''-' , ^"'^ (21) 

« + ^2 + (1^)2 sinh(6y^p2 + ^2 + (1^)2) 

where the function h{d) is given by : 

Performing a change of variables v = (p^ + + {^)^y , and now going back to the four- 
dimensional case it is possible to write hiv, a, b, d) \d=A = h{y, 0., b) as : 

cosh((fe — 2y)v) 



I °° roo 

hiy;a,b) = 7^J2 • uh 

^Txa ~^ J a smh bv 

where the lower limit of the above integral is given by 



(23) 



a = Jm? + (— )2. (24) 
V a 

It is important to keep in mind that the situation is completely different for c? 7^ 4, since after 
change of variables the term (f^ — — it will appear in the integrand of Eg. (|21|) . As 

a conseguence of this fact, it is more difficult to perform algebraic manipulations that allow us 
to analytically regularize R2{x,a,b,d) and l2{y,a,b,d) and also other expressions. Nevertheless, 



utilizing the same technigues used in refs. [^, it is possible to perform all the calculations 
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for a higher dimensional Euchdean space. This generahzation will be presented in section IV. Going 
back to Eg. (^31) and using trigonometric identities, we have : 

h{y, a, b) = hiiy, a) + I^^iv, b) + 122(2/, «, b) (25) 

where : 



hi{y,a) = ^j:exp (-2y^m2 + j , (26) 

I+{y,a,b) = - dq {cothq-l)e—, (27) 

^'Tcya Jba 



and finally 



^22(2/, b) = 1 H / dq (coth g - 1) 



222 



(2^ 



An exact expression in the massless case can be obtained from l2i{y,a). Summing the geometric 
series /2i(2/, a)|m=o can be written as 

1 1 

hi{y,a)\m=o = 2^y ■ (29) 

Airay g— - 1 

We will use the Bernoulli polynomials and numbers in the next sections. Consequently, let us 
introduce this Laurent expansion to study the behavior of l2i{y,0')\m=o in the neighborhood of 
y = 0. To find the principal part of the Laurent series of hiiy, a)|m=o around y = let us use the 
Bernoulli polynomials, which are defined by the generating function 

J. „xt 00 j.n 

^— = 5:5„(x)-, |t|<27r, (30) 
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and the Bernoulli numbers Bn = Bn{x = 0), (i?o = = ■■)■ Using this Laurent expansion 
one finds 

/2i(l/,a)Uo=(^^ + 4^ + ^E^( — ) j. (31) 

We note that we still have to calculate N{x, y, a, b). Using again the trigonometric identity given 
by Eq.(0), it is possible to write N{x, y, a, b, d) = Ni{x, a, b, d) + N2{x^ y, a, b, d) where we have : 

iVi(x, a, 6, d) = — --Wt y / ^^''"'pt " r (32) 

and 



cos 



cosh((6 - 2y)Jp^ + m2 + (H^)2) 



(33) 

sinh(6y'p2 + ^2 + (i^)2) 
Let us study the expression given by Ni{x, a, 6, d). It is possible to write Ni{x, a, b, d) = Nii{a, b, d) + 
Ni2{x, a, b, d), where A''ii(a, b, d) and Ni2{x, a, 6, d) are given by 



4(27r)'^-2a6 7 ^ (p 2 + ^2) ^ 

and 



, , . 1 [ ,d-2 1 cosh((a-2a;)^(p2 + m2)) 
iVi2(x, a, 6, rf) = - _ / d'' V , . , . ■ (35) 



The expression given by Eq.(0) can be easily calculated using dimensional regularization and can 
be written as 

^"("•'•'') = 4^K27^r(2-f)K)t- (36) 
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This contribution for the vacuum fluctuations also demands bulk conterterm. We have to deal 

with the expression of Ni2{x,y,a,b,d). Integrating over the solid angle, changing the variables 

and using the fact that we are in d = 4, i.e. defining N^i^x, a, b, d)\d=4, = Ni2{x, a, h) we have 

, ,x 1 /""^ , cosh((a - 2x)f ) , , 

Nuix,a,b) = -— dv ^. 37 

OTTO Jrri smh av 

Again, using trigonometric identities, we have : 



Ni2{x,a,b) 



±^-2xm ^ i ^^^^^ g2,f ^_ r^q (coth g - 1) e-^ 

2x 2 J am 2 J am 



■ (38) 



87ra6 

We note that Eg. (pS]) has surface divergences when x and x — > a. The structure of the 
divergences of Ni2{x, a, b) will be studied further. 

As a next step in the discussion, let us investigate the part of the one-loop vacuum fluctua- 
tions that contains edge divergences that is given by N2{x,y,a,b). Again, integrating over the 
solid angle, changing the variables and using the fact that we are in d = 4, the expression for 
N2ix, y, a, b)\d=i = N2{x, y, a, b) is given by 

1 ,2n7rx, f°° , cosh((6 - 2?/)t;) 

iV2 X, a, 6 = ^ E cos / dv \^ 1 , > 39 

4a a J a smh bv 

where the lower limit of the above integral is given by a, defined in Eg. (p4D . Using trigonometric 

identities it is possible to write N2{x, y, a, b) as 

N2 (x, y,a,b) = N21 (x, y,a) + N22 (x, y , a, 6) , (40) 

where : 

iV2i(x,y,a) = E cos(^)e-^^" (41) 
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and 



N22{x, y, g, &) = — cos( ) / (if (coth(6t>) — 1) cosh 2t>?/. (42) 

4(3 a, J a 

From now we will put m = 0. Using the Poisson's kernel given by 

oo 1 — 



l + 25]r"cosn0= -— ^, 

^ 1 — 2r cos (p + 



< r < 1, (43) 



it is posssible to obtain a closed expression for Eg. (^T]) , for points outside the boundaries. Note 
that N2i{x,y,a,b) diverges for y — > for any x. Let us study the behavior of N2i{x,y,a,b) near 
the edges x = 0,y = and also x = a,y = 0. Using the Poisson's kernel, the function N2i{x, y, a, h) 
in X = or X = a is given by 

1 1 11 

N2i{x,y,a)\^=a = — —y + — . (44) 

IQay IQay + 1 IQay + 1 

Again using the generating function of the Bernoulli numbers it is possible to find the principal 
part of the Laurent expansion of N2i{x, y, a, 6)|^=o around y = 0, and also N2i{x, y, a, b)\x=a around 
y = 0. The Laurent expansion around y = allow us to write N2i{x, y, a, &)|x=o as 

N2i{x, y, a)U=o = + ^ + ^iv)- (45) 

We note that Eg . (|4l|) has two edge divergences, one at x = y = and the other at x = a, y = 0. 
To investigate the divergences of Eg.(|42|) let us rewrite N22{x,y,a,b) as 

N22{x, y, a, b) = N+{x, y, a, b) + N:^^{x, y, a, b) (46) 
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where 



t2ix, y,a,b) = ■^£cos(^^)/ dq {cothq-1) 
oab ^ a Jab 



2gv 

e b 



(47) 



N22{x,y,a,b) 



Sab 



1 



cos( ) / aq [coth q — I) e 



(48) 



In the next section, we will show that N22{x,y,a,b) has surface divergences when |/ — > 6 for all 
X, and also two edge divergences, one at x — 0, y — b and the other at x — a, y — b. We note 
also that N22{x, y, a, b) is finite everywhere. Using a combination of dimensional and zeta function 
regularization, in the next section we will present the general method to regularize the one-loop 
vacuum field fiuctuations in the rectangular waveguide. 

3 Analysis of the surface divergences in the one-loop vac- 
uum fluctuations. 

The purpose of this section is to present a general method to analytic regularize the one-loop 
vacuum fiuctuations associated with the confined scalar field. We first present the structure of the 
surface and edge divergences of the vacuum fiuctuations associated with a massless scalar field in 
a four dimensional Euclidean space. As we discussed in the previous section, it is possible to write 
the one-loop vacuum fiuctuations as 



The first expression that we have to deal with is T{a,b,d). As we discussed in the previous 



Tdd{x, y, a, b, d) = T{a, b, d) + T{x, a, b, d) + T{y, a, b, d) + T{x, y, a, b, d) 



(49) 
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section, the analytic structure of T(a, 6, d) was carefully analysed by Kirsten |Q and it is not 



necessary to repeat the calculation again. The second term of the above equation and third one 
can be written respectively in terms of Ri{a, b, d), R2{x, a, b, d), Ii{a, b, d) and finally l2{y, a, b, d). 
The polar structure of Ri{a,b,d) and Ii{a,b,d) can be found in the literature, and we will not 
repeat the analysis that was done in these papers. The next quantity that we have to regularize is 
R2{x, a, b, d) and also hiy, a, b, d). Since both cases are equivalent, let us study only the expression 
given by R2{x, a, b, d). It is instructive to study first the simpler case b ^ a which will give us an 
indication of the behavior that may be encounter. Let us now proceed with the calculations in a 
general d-dimensional Euclidean space. Defining R2{x, a, b, (i)|b>a = r2(x, a, d), we have : 



1 f 1 cos}i{{a — 2x) \l fP' + m?) 



r2(x,a,d) = ;j—r / d'^ ^p—=== , ^ . (50) 



We will again use the fact that d p = p dpdVtd-i and / dVtd-i = f^iEiy- Let us choose now 



m = 0. Defining h2{d) by 



h2(d) = , (51) 

21 J 2^-%^r(^)' ^ ' 



it is possible to write r2(x, a, d)\m=o as 



r2{x,a,d)\m=o = ^h2{d) 

oo 



PCX) 

/ dk A;'^"^(coth ka - 1) cosh 2kx 
Jo 



+ / dk /t'^-^g-a/cx 







(52) 



In the first integral for large k, (cothfca — 1) has the behavior: (cothfca — 1) ~ e Moreover, 
the second integral in the above equation is ultraviolet finite for a; 7^ 0. Let us define t = ka and 
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q = kx in the first and second integrals above respectively. Then Eq.(^) becomes 



r2{x,a,d)\ra=o 



2a<^ 
1 



1 Ixt 
— r/i2(rf) / dtf^-^icoiht - 1) cosh( ) 



(53) 



The second term in the above equation gives us the well known result that for a massless scalar 
field {ip^i^x)) diverges as (in a four- dimensional space) as we approach the plate 0. In order to 
analyze the behavior of r2(x, a, d) \m=o around x = a, let us use the following integral representation 
of the gamma function, 



poo 2 

dtf-^e-''* = —T{^j), Re{fi)>0, Re{u) > 



(54) 



and also the following integral representation of the product of the Gamma function times the 
Hurwitz zeta function 



/ c?tt^-^e-"*(cotht - 1) = 2^-^r(/i)C(;U, - + 1) Re{a) > 0, Re{fi) > 1, 
Jo 2 



where C{s,u) is the Hurwitz zeta function defined by |64 



(55) 



C(^,«) = E 



+ 



Re{s)>l, 0,-1, -2. 



Then, using Eqs. (M) , (p5|) and (|56D in Eq.(^) we have that : 



r2{x,a,d)\^=o 



2 a 
+ 7:r4^/i2(rf)r(rf-2) 



X 



X 



2^-^r(rf - 2) C(rf - 2, - + 1) + Cid - 2, — + 1) 

a a 



(2x; 



(56) 



(57) 
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Using the definition of tlie Hurwitz zeta function, it is evident that : 

^ ^c(rf-2,- + i) + c(rf-2,-- + i; 



1 ^ 1 1 1 ^ 1 , 



n=0 (n + (1 + f (« - ^)'-' n^l + (1 - f ))' 

We see that the regularized r2(x, a, (i)|m=o has two poles of order (c? — 2) in x = and in x = a. 
Note that the residues of the poles in a; = and in x = a are a- independent. The same analysis 
can be done for 12(1/, a, d) assuming a ^ b. In the next section the general case will be studied. 
Let us finally analyze N2i{x, y, a) and N22{x, y, a, b), given respectively by 

^2i(a;,2/,a) = -— 2^ cos( )e - (59) 

Say a 

and 

^22(3:, 2/, a, 6) = — — V cos( )/ (cothg - 1) cosh(— — ). (60) 

4ao 1 a J^b b 

To find the analytic structure of N2i{x, y, a) we can expand the general term in the sum in power 
series, commute the two summations, and use analytic continuation in the zeta function that will 
appear. The process will in general produce an extra term, which is generated by commuting the 
convergent exponential summation J2m with the new divergent summation J2k (fo^' details see for 



e.g. refs. [0 |^ HHI). In our case this term vanishes due to the power of n. Let us express the 



sum that appears in Eq.(^) in terms of the complex variable z = ix — y: 



1 ,2mTZ,] , , 

iV2i X, y, a) = -^Re ^ exp . 61 

Say a J 
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Expanding around z = 0, (using the Bernoulli expansion) will produce : 

We see that the edge divergence appears in the term 1/z. Taking the real part : 

N2^{x,y,a) = - I - ^Re{f,{z)} , (63) 

where fi{z) is an entire function of z and is given by : 

Expanding around z = ia will produce : 

;V„(..,.a) ^ -2.ne [t (!^)' (i^C(-*) - (^) (. - «)-'). (65) 

Taking the real part of 1/(2; — ia), we have : 

N2i{x,y,a) = -——^^ r^-^f2{z), (66) 

loTT [y^ + [x — a)^) Say 

where f2{z) is also an entire function of z, and is given by: 

/.W^-{g(a.)'(i^CM.,}. 

To find the analytic structure of N22{x, y, a, b), it is enough to analyse the quantity N22{x, y, a, b), 
which is given by the Eq. (^7D (the other quantity A^22(^5 ^) is finite everywhere). To calculate 
the integral in Eq.(^7D, we can express coth(6f ) in exponential functions and expand the integrand 
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in power series. The integral can be easily evaluated, and the result is given by 

N+(x,y,a,b) = r ^ cos ^. 68 

This result can be written in terms of the Lerch's trancendent function $(2;, c, v) which is defined 
as follows [[7^ : 

which is valid for \z\ < 1. By analytic continuation it can be extended to the whole complex plane. 
This function has singularities at z = 1 and c = or c = 1 and when f is a non-positive integer 
and Re(c) is also non-positive. We note that only the term A; = has a surface divergence at 
y = b. The remaining part of the series in k is finite everywhere and we call it F{x, y, a, b). Using 
the same procedure used to analyse N2i{x, y, a), we can define z = {b — y) + ix and expand the 
term A; = around z = or around z = ia. Therefore we have 

N22ix, y, a, b) = B{x, y, a, b) + F{x, y, a, 6), (70) 

where the first expansion B{x, y, a, b) is given by 

^' ^) = - 16n{x^Hy-by) - M^^^^'^' ^^'^ 
where hi{z) is an entire function given by 
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We note that in this case B{x,y,a,b) has a edge divergence for x = 0, y = b. For the second 
expansion, B{x, y, a, b) is given by 



We note that in this case B{x, y, a, b) has a edge divergence for x = a, y = b. A general picture 
that emerges from the above discussion is the following: we have found that in order to eliminate 
the ultraviolet divergences of the one-loop vacuum fluctuations we have to introduce not only 
the usual bulk counterterms, but also counterterms concentrated on the boundaries. In the next 
section we will show how it is possible to generalize our local analysis for a massive field defined 
in a higher dimensional Euclidean space. 



4 The one- loop vacuum fluctuation in the waveguide in 
a higher dimensional Euclidean space. 



In this section we will discuss how to generalize our calculations to a higher dimensional Euclidean 
space. We are interested in investigating the one-loop vacuum fluctuations associated with a 
massive scalar field defined in a generic d-dimensional Euclidean space, where d — 2 coordinates 
are unbounded while the first two lie in some finite interval. As discussed in the previous sections. 




(73) 



where h2{z) is an entire function given by 




(74) 
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we found that the vacuum fluctuations can be expressed as : 

(y{x,y)) = ^T{a,b,d) + ^{Ri{a,b,d) + Ii{a,b,d)) + 

i {R2{x, a, b, d) + hiy, a, b, d)) + N{x, y, a, b, d). (75) 

In the previous sections, we studied the expression given by T(a, 6, rf), i?i(a, 6, d) and also /i(a, b, d) 
in the massless case, which demands bulk counterterms. In the following, we will present a detailed 
analysis of the contributions of the vacuum fluctuations that demand counterterms concentrated 
on the boundaries. We will show that some of the terms of the one-loop vacuum fluctuations, can 
be expressed in terms of Bessel functions and also Hurwitz zeta functions. 

From the previous analysis for the massless field in a four-dimensional Euclidean space, we 
found that to perform our calculations in a generic d-dimensional space, for the massive field, it 
is natural to write /2(?/, a, d) in the following way : 

h {y, a,b,d) = I21 (y, a, d) + I22 {y, a, b, d) , (76) 

where the functions /2i(y, a, d) and 122(1/, b, d) are given by : 

l2dy,a,d) = -h{d)Y. dv{v^-m'-{—f)^e-'^y, (77) 

and 

1 °° roo jlTl d-4 

l22{y,a,b,d) = -h{d)Y, dv{v^ - - { — )2) — (cothfew - 1) cosh2w2/. (78) 

a .^^^ J a a 

The lower limit of both integrals is given by a defined in Eg . (|2^) and the coefficient h{d) is an 
entire function defined in Eg. ([2^) . Using the following integral representation of the modified 
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Bessel functions of third kind, or the Macdonald's functions K^{x) : 

/ {x^ -u'Y-^e-^''dx = ^{—y-^T{v)K,^i{un), (79) 

which is vahd for for m > 0, -Re /x > and Reu > 0, hiiv, a, b, d) can be written in terms of these 
functions. A simple substitution gives 

h,{y,aAd) = -77r^^T.i-)'^K^{2ay). (80) 

« i^VT^r t^i y 

Equation (|80|) can not be evaluated exactly, even for the massless case, but for small arguments of 
the Bessel function it is possible to show how does l2i{y, a, d) contributes to the surface divergences 
of the one-loop vacuum fluctuations. 

The expression 122(1/; '^5 b, d) can also be computed in the both cases, the massive and the mass- 
less one. To ilustrate the method that we discussed in previous sections, let us present the massless 
case calculation. The massive case follows the same procedure. To calculate 122(1/, b, d)\m=o note 
that after change of variables 122(2/, a, (i)|m=05 contains a power of a binomial. When d is even 
(for d > 4) the power is an integer and the use of the Newton's binomial theorem will give a very 
direct way to generalize our further results. When d is odd {d > A), the expansion on the binomial 
yields an infinite power series. The same technique was used in refs. It is worthwhile 



to remark that the study of even case follows the same procedure, only with a finite numbers of 
terms in the binomial expansion. Consequently, let us study the most interesting case, i.e the odd 
dimensional case. Let us use the following series representation for (1 — (— )" 



d-4 

■ va ' 



l~{-fY-^ = ±{-lfCiC^?\ (81) 



va ' ' ^ ~ ' ^ va 



27 



where the are the generahzations of the binomial coeficients, given by: (7° = 1, = ^, = 
^p{p — I),.-- Cp = yipip ~ ~ k + \)) and p = Note that the generahzation of the 

binomial series is valid even for any complex exponent p. In other words, for f > ^ we have 
an everywhere convergent power series in p, hence a continuous function on p in the complex 
d plane. In the following, it is convenient to define the quantities: C^^\d,k) = {—l)^Cph{d), 
C^^^{d, k) = 7c^''C^^\d, k) and finally C(3)(d, k) = ^''^^Ztl^^ C^^^d, k). After a change of variables, 
using Eq.(^ in Eq.(|7S]), /22(y; a, c?)|m=o becomes 

h2{y, a, 6, d)U=, = -— Y: C^'\d, k){-f' L u'-'-^^coihu - 1) cosh(^), (82) 

where C^'^\d, k) is an entire function in the complex d plane. The integral that appear in Eq. 
cannot be evaluated explicity in terms of well known functions. Nevertheless it is possible to write 
Eq.([8^) in a convenient way where the structure of the divergences near the plate when y ^ b 
appear. Let us split the 122(2/5 0,, b, d)\m=o in the following way : 

/22(y, a, 6, d)\m=o = l22iy, 6, d) + I^^iy^ b, d), (83) 

where 

I<(y,a,b,d) = -— Y: C7(2)(d,A;)(-)2'=5:n2M^^^rf««'^-^^2'=(coth«-l)cosh(-^), (84) 

k=0 ^ n=l ''— ^ 

and 

I>(y,a,b,d) = —^ Y C^'\d,k)i-rY^"' J^duu'^-''-'''icothu-l)cosH^). (85) 
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Although the integral that appears in Eq.(^) can not be expressed in terns of known functions, the 



leading divergences of l22iy, cl, b, d)\m=o are contained in the Eq. (|8^) . Consequently, let us analyze 
the divergences that appear in Eq. (^) . Using the fact that f{u)du = f{u)du — Jq f{u)du, 
we reexpress Eq.(P3) as 



I^^iy, a, b, d) = 122(2/, a, 6, d) + I^^iy, a, 6, d), (86) 

where 

,2uy 



7) 



1^2 h °o poo ' 

nJy,a,b,d) = —=r^ V C^^\d,k){-f^y n''' duu''-^-^\coi\iu - I) cosh(- 
and 

1 h °° r— 91111 

I'^,iy,a,b,d) = --^ y C^^\d,k){-y'yn^'' " duu''-''-^\cothu-l)cosH^). ^ooj 
a¥ 3 a Jo b 

Expanding cosh{^) in powers series, the integral that appears in /|2(2/, b, d) is a Debye integral. 
Although this can be evaluated, the study of /f2(2/, a, c?) is sufficient to obtain the divergences 
of /^(y, a, 6, (i). Again, we note that the integral that appear in Eq . (pT]) can be written in terns 
of products of zeta. Gamma function and Hurwitz zeta function. Thus we have : 

12 00 :121b ^ 



I<{y,a,b,d)\m=o = -^-TT E C^^\d,k){-f^yn^^ ( ' duu^'^'^^cotliu - I) coshC- 



b 



l;^C^'\d, 0) (ad - 3, -I + 1) + ({d - 3, 1 + 1) V (89) 



Since C^^\d,0) is an entire function in the complex d plane, it is clear that we have the same 
surface divergences that we studied before. To study the massive case 122(2/, a, c?) we have to 
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use that same procedures. Let us show how our procedure can be systematically applied in the 
other expressions. From the previous calculations we express the contribution to the vacuum 
fluctuations that contains edge divergences giving by N{x, y, a, b, d) in the following way : 

A^(x, a, 6, d) = Nii{a, b, d) + Ni2{x, a, 6, d) + N2{x, y, a, b, d), (90) 

where the Nii{a,b,d) term contains contributions to the vacuum fluctuations that demands bulk 
counterterms in order to renormalize the one-loop vacuum fluctuations. The next contribution 
that we have to study is Ni2{x, a, b, d). Starting from Eg . (|35D integrating over the solid angle and 
changing variables, it is also convenient to split Ni2{x, a, b, d) in the following way : 

iVi2(x, a, 6, d) = N^2ix, b, d) + Nl^{x, a, 6, d), (91) 

where each of the terms of Eg . (^TD are given respectively by 

Nl^{x, b, d) = --Th{d) / dv {v^ - m^)^e-^''\ (92) 



and 



1 r°° d-A 

Nl2{x,a,b,d) = -—h{d) dv {v'^ - m^)~ {cothav - I) cosh.2vx. (93) 

40 Jm 



We can again use the integral representation of the modifled Bessel functions of the third kind 
given by Eg. ([79|) , to write Eg.(p2D as 

NUxAd) = -i^^-^ip3^(^)^ir^(2mx). (94) 
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To calculate N^2i^y d)-, let us study first the massless case. In this case, first changing variables 
and using Eq.(^5|) we get 



NUx, a, 6, d)U=o = -^^r(rf - 3) (Cid - 3, - + 1) + C(rf - 3, -- + 1)) . (95) 

Next we consider the massive case. Note that we have exactly the same situation discussed before. 
After a change of variables, Ni2{x,a,b,d) also contains a power of a binomial. Let us discuss 
again the most interesting situation, the odd dimensional case. First, let us the same power series 

2 d — 4 

expansion that we used before for the quantity (1 — given by 

2 oo 

rr^^^ _ (96) 
Substituting Eq.(^) in Eq.(|93D and again changing variable, we have 

k=0 

A natural way to obtain N'l2{x,a,b,d) expressed in terms of the Hurwitz zeta function is the 
following. Let us split N^2i^y d) in the following way : 



N'^2ix,a,b,d) = -—l-^h{d)Y.{-lfCl{maf''j du u''-^-^\coth u - I) cosh (97) 



NI2 (x, a, 6, d) = (x, a,b,d) + (x, a, 6, d) , (98) 

where 

N^<{x,a,b,d) = -^-^ J2 C^^\d,k){amf'' J du u''-^-^\coth u - 1) cosh{-^) , (99) 



k=0 

and 



1 °o roo lnT 

N'>{x,a,b,d) = -^-^ ^ C^^\d,k){amf'' J ^^^'^-^-^'^(cothu - 1) cosh(-^). (100) 
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We have the same situation that we studied before. Let us investigate first the N\^{x,a,h,d). 
Using that f{u)du = f{u)du — Jq f{u)du, it is possible to write Eq.(^) in the following 
way: 



Nl<{x,a,h,d) = -—— Y: C^'\d,k)i—y''(cid-3-2k, — + l) + ad-3-2k,- + l) 



1 2 i-ain Out 

y C'^^\d,k)(am?'' duu''-^-''-(cothu-l)cosh(^), (101) 
zrt, JO a 



k=0 

where the singularities of N^2{^^ ^i ^i appear at x ^ a. We have to analyze the A''j'^(x, a, 6, d). It 
is clear that in a even dimensional Euclidean space for x < a the integral that appear in Eq. (|100|) 
is convergent. The odd dimensional case also can be studied. In this case we have to expand 
cosh{^^) in power series and it is clear that in the new integral that appears in Eq.( |100| ) demands 
a generalization of the Debye integrals. The calculations that we presented can also be used to 
write N2{x, y, a, b, d) in the same form that we obtained for the other expressions. Unfortunatelly 
one cannot perform the summation in n. For completeness, we will only write the remaining 
expressions. Let us start from : 

N2{x,y,a,b) = N2i{x,y,a,d) + N22{x,y,a,b,d), (102) 

where each term of the above equation is given by 

1 , / ,x ,2mxx, f°° , , o o ,n7r d-i 

n=l 

and 



N2i{x, y, a, d) = —h{d) y cos(^) T dv{v' - - {—f)^e-^^y, (103) 
27ra rr^ a Ja a 



1 °° 2mTX 

N22{x,y,a,b,d) = - — /i(rf) ^ cos(- 



2™ a 
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/ dv{v'^ - - { — Y)~ (coth bv - I) cosh 2vy. (104) 

Ja a 

Note that the lower hmit of both integrals are given by a defined in Eg. (P^ and h{d) was defined 
in Eq.(^). Let us present the N2i{x,y,a,b,d) in the massless case. Using the same integral 
representation of the Bessel functions that we used before gives 

d-3 

2 ^ 2mTX ( mx\ ^ 2mTy 
N21 (x, y, a, d) \m=o = ,^ , , , 2^ cos( ) — Kd^[ ). (105) 

The massive case can also be presented but in both cases it is not possible to perform the sum- 
mation in n. Finally, let us present the expression for N22{x, y, a, b, d) in the massless case and for 
the odd dimensional case : 

N22ix,y,a,b,d)\m=o = ^^"^ E C^'H^, fc)!-^ E cos(^) 

duu''~^-^\cot\iu - 1) cosh(^). (106) 

a 

The same procedure that we used in l22{y,a,b,d), can be repeated again. Unfortunatelly one 
can not perform the summation in n that appears in Eq.( |106| ). Nevertheless the same analysis 
that we did in the end of section III can be performed to study the edge divergences that appear 
in the one-loop vacuum fiuctuations associated with the scalar field confined in the interior of a 
waveguide in a higher dimesional Euclidean space. 
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5 Discussions and conclusions 

In this paper, we first obtained the regularized one-loop vacuum fiuctuations associated with a 
massless scalar field defined in the interior of an infinity waveguide of rectangular cross section in 
a four-dimensional Euclidean space. Then, we discussed how it is possible to generalize our results 
for a massive field defined in a higher dimensional Euclidean space. 

Let us summarize our motivations and the results. First, in the rectangular waveguide con- 
figuration the calculations can be done analytically, and since in this geometric configuration the 
electromagnetic field can be described by two massless scalar fields with Dirichlet and Neumann 
boundary conditions, our calculations can be used to describe the vacuum fiuctuations associated 
with the electromagnetic field in the interior of the waveguide. Also, the rectangular geometries 
are very convenient to carry out experiments on the effect of boundaries on the atomic levels of 
atoms and on the rate of spontaneous emission. It is well known that in waveguides, it is pos- 
sible to obtain situations where the spontaneous emission of atoms can be suppresed and also 
enhanced. Finally, the calculation of the regularized one-loop vacuum fluctuations in the inflnite 
waveguide of rectangular cross section using the combination of dimensional and zeta function 
analytic regularization has not been discussed in the literature, at least as far as we know. 

We first rederive a well-know result that surface and edge divergences appear in the one-loop 
vacuum fluctuations as consequence of the uncertaintly principle. There are at least two different 
possible solutions that can eliminate these divergences. The first is to take into account that 
real materials have imperfect conductivity at high frequencies. As was stressed by many authors, 
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the infinities that appear in renormahzed values of local observables for the ideal conductor (or 
perfect mirror) represent a breakdown of the perfect-conductor approximation. A wavelength 
cutoff corresponding to the finite plasma frequency must be included. The second would be 
given by a quantum mechanical treatment of the boundary conditions It was shown |53 

that position fluctuations of a reflecting boundary also remove divergences in the renormalized 
values of local observables at least in the fiat plate configuration. A logical going is to use an 
analytic regularization procedure to identify these divergent terms. As we discussed counterterms 
concentrated on the boundaries produce a finite one-loop vacuum fluctuations in the interior of 
the waveguide. A natural extension of this paper is to go beyond the one-loop approximation, 
investigating interacting field models in the presence of a waveguide in a d-dimensional Euclidean 
space. This topic is under investigation by the authors. 
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A The one-loop vacuum fluctuations and the generahzed 
zeta function method 

In aim of this appendix is to discuss the hnk between the one-loop vacuum fluctuations and the 



generalized zeta function method |3J] as a standard formalism to regularize products and 
determinants. Since our prime goal in section IV was to calculate the one-loop vacuum fluctuations 
in the massive model, we will present the formal relation between the one-loop vacuum fluctuations 
and the effective action in the one-loop approximation. We are following the treatment used by 
Dittrich and Renter [0]. 

Let us consider the generating functional of the complete Schwinger functions for a scalar field 
in a d-dimensional Euchdean space: 

Z{J) = I /}^e-^[^l+/'^'^-^(^)^(^), (A.l) 



where Dip is the appropriate measure, and S[ip\ is the classical action associated with the scalar 
field. The quantity Z{J) can be regarded as the functional integral representation for the imagi- 
nary time evolution operator {(p2\U {t2,ti)\ipi) with the boundary conditions ip{ti,x) = fi{x) and 
(f{t2, x) = (p2{x). The quantity Z{J) gives the transition amplitude from the initial state \ipi > to 
a final state |(y92 > in the presence of some scalar source of compact suport. As usual W{J), the 
generating functional of the connected correlation functions shall be given by W{J) = lnZ{J). In 
a free theory the partition function Z{J) and also W{J) can be calculated exactly. In the presence 
of the waveguide we must assume that the path integral must be taken over the space of functions 

36 



that vanish in the boundaries of the waveguide. Since our aim is to discus the non-interacting 
theory, let us assume that J — 0. Thus, defining the d-dimensional laplacian by A, we have: 

J Dirichlet 

"Dirichlet" in the path integral means that we are performing the path integral over functions 
that vanish on the boundaries. Thus the generating functional can be written as 

Z( J)| j=o = N{det{-^ + m2))-5, (A.3) 

where A?" is a normalization factor that does not contributes to the free energy. For simplicity, let 
us discuss the four-dimensional case. The global generahzed zeta function is defined by 

C-A-,m<s) ^ N' J dkJ dh E {kl + kl + i—r + A' + ruY', (A.4) 

n,n'=l ^ 

where N' is a normalization constant. It follows that 

Z(J^O) = expC-CA+^m, (A.5) 
where C'(0) = ;^C('S)|s=o and the operator (—A -|- m^) has the spectrum 

{kl + kl + i—f + &^ + m^), ko, ki e ^, n, n' e N. (A.6) 
a 

We can use the same procedure to define the local generalized zeta function which is related with 
the effective action. The local generahzed zeta function is given by 
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where x and y are two cartesian coordinates of one point of the infinite waveguide with rectangular 
crosss section. Thus we have the formal expression: 

{-m^r{^'{x,y)) = as,x,y). (A.8) 

To give a precise meaning of the derivative one can use for example the Liouville's concept of 



fractional derivative |75|. This result show that the one- loop vacuum fluctuations determines the 



effective action in the one-loop level. 
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